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Abstract. 

Recently, a new non-Standard Model neutrino interaction mediated by a light scalar 
field was proposed, which renders the big-bang relic neutrinos of the cosmic neutrino 
background a natural dark energy candidate, the so-called Neutrino Dark Energy. 
As a further consequence of this interaction, the neutrino masses become functions 
of the neutrino energy densities and are thus promoted to dynamical, time/redshift 
dependent quantities. Such a possible neutrino mass variation introduces a redshift 
dependence into the resonance energies associated with the annihilation of extremely 
high-energy cosmic neutrinos on relic anti-neutrinos and vice versa into Z-bosons. In 
general, this annihilation process is expected to lead to sizeable absorption dips in 
the spectra to be observed on earth by neutrino observatories operating in the relevant 
energy region above 10 13 GeV. In our analysis, we contrast the characteristic absorption 
features produced by constant and varying neutrino masses, including all thermal 
background effects caused by the relic neutrino motion. We firstly consider neutrinos 
from astrophysical sources and secondly neutrinos originating from the decomposition 
of topological defects using the appropriate fragmentation functions. On the one hand, 
independent of the nature of neutrino masses, our results illustrate the discovery 
potential for the cosmic neutrino background by means of relic neutrino absorption 
spectroscopy. On the other hand, they allow to estimate the prospects for testing its 
possible interpretation as source of Neutrino Dark Energy within the next decade by 
the neutrino observatories ANITA and LOFAR. 



Probing Neutrino Dark Energy with Extremely High-Energy Cosmic Neutrinos 2 
1. Introduction 

According to Big Bang Cosmology, in an expanding universe the freeze-out of a 
particle species occurs, when its interaction strength is too small to keep it in thermal 
equilibrium. Neutrinos, being the particles with the weakest known interactions, 
therefore are assumed to have already decoupled when the universe was just ~ 1 s 
old, thereby guaranteeing a substantial relic neutrino abundance today with average 
number density n Voi = n Poi = 56cm~ 3 per neutrino species i = 1,2,3. Yet, in turn, 
the weakness of the neutrino interactions so far has spoilt all attempts to probe this 
1.95 K cosmic neutrino background (Cz/B), which is the analog of the 2.73 K cosmic 
microwave background (CMB) of photons, in a laboratory setting [1-4]. However, other 
cosmological measurements, such as the light element abundance, large scale structure 
(LSS) and the CMB anisotropies are sensitive to the presence of the Cz/B and therefore 
have provided us at least with indirect evidence for its existence (see e.g. Ref. [5] for a 
review) . 

Independently, Type la Supernova (SNIa) results (e.g. [6]), supported by CMB [7] 
and LSS data (e.g. Refs. [8,9]), strongly suggest the existence of an exotic, smooth 
energy component with negative pressure, known as dark energy, which drives the 
apparent accelerated expansion of our universe. Recently, Fardon, Nelson and 
Weiner [10] have shown that the relic neutrinos, which constitute the Cz/B, are promoted 
to a natural dark energy candidate if they interact through a new non-Standard-Model 
scalar force - an idea which has great appeal. Neutrinos are the only Standard Model 
(SM) fermions without right-handed partners. Provided lepton number is violated, 
the active (left-handed) neutrinos are generally assumed to mix with a dark right- 
handed neutrino via the well-known seesaw mechanism [11-14], thus opening a window 
to the dark sector. Therefore, it would not seem to be surprising if neutrinos, whose 
interactions and properties we know comparably little about, were sensitive to further 
forces mediated by dark particles. Moreover, the scale relevant for neutrino mass squared 
differences as determined from neutrino oscillation experiments, 5m u 2 ~ (10 -2 eV) 2 , is 
of the order of the tiny scale associated with the dark energy, (2 x 10 -3 eV) 4 . 

As a consequence of the new interaction in such a scenario, an intricate interplay 
links the dynamics of the relic neutrinos and the mediator of the dark force, a light scalar 
field called the acceleron. On the one hand, the neutrino masses m u . are generated 
by the vacuum expectation value A of the acceleron, m Vi {A). Correspondingly, the 
A dependence of the masses m u .(A) is transmitted to the neutrino energy densities 
PvA m Vi(A)) since these are functions of m Vi (A). On the other hand, as a direct 
consequence, the neutrino energy densities p Vi (m Vi) A) can stabilize the acceleron by 
contributing to its effective potential V c s(A, p Ui ), which represents the total energy 
density of the coupled system. Moreover, cosmic expansion manifests itself in the 
dilution of the neutrino energy densities p Vi (z) ~ (1 + z) 3 . Therefore, it crucially 
affects the effective acceleron potential V^(A, p Vi (m Vi , z)) by introducing a dependence 
on cosmic time, here parameterized in terms of the cosmic redshift z. For a homogeneous 
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configuration, the equilibrium value of the acceleron instantaneously minimizes its 
effective potential V e g(A, p Ui (m Vt , z)) and therefore also A(z) varies on cosmological time 
scales. Finally, in turn, since the neutrino masses m Vi {A) are sensitive to changes in A, 
they are promoted to dynamical quantities depending on z, m u .(z), that is depending 
on cosmic time. To summarize, the variation of the neutrino masses represents a clear 
signature of the so-called Neutrino Dark Energy scenario. 

In a subsequent work Fardon, Nelson and Weiner [15] presented a supersymmetric 
Mass Varying Neutrino (MaVaN) model in which the origin of dark energy was 
attributed to the lightest neutrino v 1 and the size of the dark energy could be expressed 
in terms of neutrino mass parameters. By naturalness arguments the authors concluded 
that the lightest neutrino still has to be relativistic today, thereby preventing potential 
instabilities [16-18] which could occur in highly non-relativistic theories of Neutrino 
Dark Energy. 

The rich phenomenology of the MaVaN scenario has been explored by many authors. 
The cosmological effects of varying neutrino masses have been studied in Refs. [19,20] 
and were elaborated in the context of gamma ray bursts [21]. Apart from the time 
variation, the conjectured new scalar forces between neutrinos as well as the additional 
possibility of small acceleron couplings to matter lead to an environment dependence of 
the neutrino masses governed by the local neutrino and matter density [10,22,23]. The 
consequences for neutrino oscillations in general were exploited in Refs. [22,24] and in 
particular in the sun [25-27], in reactor experiments [26,28] as well as in long-baseline 
experiments [29]. 

In light of the possible realization of Neutrino Dark Energy in nature, a (more) 
direct detection of the Cz/B should be thoroughly explored with special attention turned 
to possible new physics beyond the SM. By this means, a time evolution of neutrino 
masses could be revealed which would serve as a test of Neutrino Dark Energy. In 
addition, the general importance of a (more) direct evidence for the existence of the 
Cz/B lies in a confirmation of standard cosmology back to the freeze-out of the weak 
interactions and therefore thirteen orders of magnitudes before the time when photons 
where imprinted on the last scattering surface. 

An appealing opportunity to catch a glimpse of the Cz/B as it is today emerges 
from the possible existence of extremely high-energy cosmic neutrinos (EHECz/'s). Such 
EHECz/'s can annihilate with relic anti-neutrinos (and vice versa) into Z bosons, if their 
energies coincide with the respective resonance energies E™f of the corresponding process 
vv — > Z [30-37]. These energies, 



in the rest system of the relic neutrinos, are entirely determined by the Z boson mass 
Mz as well as the respective neutrino masses m Vi . An exceptional loss of transparency of 
the Cz/B for cosmic neutrinos results from the fact that the corresponding annihilation 
cross-section on resonance is enhanced by several orders of magnitude with respect to 
non-resonant scattering. As a consequence, the diffuse EHECz/ flux arriving at earth is 




(1) 
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expected to exhibit absorption dips whose locations in the spectrum are determined by 
the respective resonance energies of the annihilation processes. Provided that the dips 
can be resolved on earth, they could produce the most direct evidence for the existence 
of the Cz/B so far. Furthermore, as indicated by the resonance energies in Eq. (jTJ, 
the absorption features depend on the magnitude of the neutrino masses and could 
therefore reflect their possible variation with time. Moreover, they are sensitive to the 
flavor composition of the neutrino mass eigenstates as well as to various cosmological 
parameters. Accordingly, the possibility opens up to perform relic neutrino absorption 
spectroscopy as an independent means to probe neutrino physics and cosmology. 

The existence of EHECz/'s is theoretically well motivated and is substantiated by 
numerous works on possible EHECz/ sources of astrophysical nature (bottom-up) (see 
e.g. [38] for a review) or so-called top-down sources (see e.g. Ref. [39] for a review). In the 
latter case, EHECz/'s with energies well above 10 11 GeV are assumed to be produced in 
the decomposition of topological defects (TD's) which originate from symmetry breaking 
phase transitions in the very early universe. 

Depending on the underlying EHECz/ sources the EHECz/ fluxes could be close 
to the current observational bounds set by existing EHECz/ observatories such as 
AMANDA [40] (see also Ref. [41]), ANITA-lite [42], BAIKAL [55], FORTE [43], 
GLUE [44] and RICE [45] which cover an energy range of 10 7 GeV < E < 10 17 GeV 
(cf. Fig. Q). Promisingly, the sensitivity in this energy range will be improved by 
orders of magnitude (cf. Fig.Q) by larger EHECz/ detectors such as ANITA, EUSO [46], 
IceCube [47], LOFAR [48], OWL [49], SalSA [51] and WRST [48] which are planned 
to start operating within the next decade (cf. Fig. [TJ). Accordingly, the prospects of 
confirming the existence of the Cz/B by tracking its interaction with EHECz/'s have 
substantially improved since the original proposal in 1982 [30]. Moreover, in the likely 
case of appreciable event samples the valuable information encoded in the absorption 
features of the EHECz/ spectra could be revealed within the next decade (cf. Fig. Q), 
rendering the theoretical exploration of relic neutrino absorption spectroscopy a timely 
enterprise. 

Note that the scenario introduced above has also attracted attention for another 
reason than the possible detection of the Cz/B- namely for the controversial possibility 
of solving the so-called GZK-puzzle to be discussed briefly in the following. Beyond 
the Greisen-Zatsepin-Kuzmin (GZK) energy, Eqzk = 4 x 10 10 GeV, ultra-high energy 
nucleons rapidly lose energy due to the effective interaction with CMB photons 
(predominantly through resonant photo-pion production) [56,57]. In the so-called Z- 
burst scenario, the observed mysterious cosmic rays above Eqzk were associated with 
the secondary cosmic ray particles produced in the decays of Z bosons. The latter were 
assumed to originate from the neutrino annihilation process outlined above [58-63]. 

However, recently, ANITA-lite [42] appears to have entirely excluded the Z-burst 
explanation for the GZK-puzzle at a level required to account for the observed fluxes of 
the highest energy cosmic rays. We would like to stress, that this only means that the 
GZK-puzzle stays unsolved (if there is any) . Moreover, this neither restricts the possible 
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Figure 1. Current status and next decade prospects for EHEC^ physics, expressed 
in terms of diffuse neutrino fluxes per flavor, F — F Va + F Pa , a — e, /j,,t. The 
upper limits from AMANDA [40], see also Ref. [41], ANITA-lite [42], FORTE [43], 
GLUE [44], and RICE [45] are plotted. Also shown are projected sensitivities of 
ANITA [42], EUSO [46], IceCube [47], LOFAR [48], OWL [49], the Pierre Auger 
Observatory in v e , modes and in v T mode (bottom swath) [50], SalSA [51], and 
WSRT [48], corresponding to one event per energy decade and indicated duration. 
Also shown are predictions from astrophysical Cosmic Ray (CR) sources [52], from 
inelastic interactions of CR's with the cosmic microwave background (CMB) photons 
(cosmogenic neutrinos) [52,53], and from topological defects [54]. 



success of producing evidence for the Cz/B by means of detecting absorption dips in the 
EHECz/ spectra nor does it spoil the possibility of gaining valuable information from 
performing relic neutrino absorption spectroscopy. 

The goal of this paper is to carefully work out the characteristic differences in 
the EHECz/ absorption features which result from treating the neutrino masses as time 
varying dynamical quantities in comparison to constants. In our analysis, we incorporate 
the full thermal background effects on the absorption process whose impact grows for 
smaller neutrino masses [36,37]. This means, that in general relic neutrinos cannot 
be assumed to be at rest. Instead, they have to be treated as moving targets with a 
momentum distribution, if their mean momenta turn out to be of the order of the relic 
neutrino masses. 

We illustrate our results for the diffuse neutrino fluxes to be observed at earth 
firstly by considering astrophysical EHECz/ sources. Secondly, we calculate the 
neutrino spectrum (both for varying and constant neutrino masses) expected from the 
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decomposition of a topological defect using the appropriate fragmentation functions 
and including the full thermal background effects. By this means, for a given EHECz/ 
injection spectrum, we present all technical tools to interpret EHECz/ absorption dips 
as soon as they are observed at earth. Thereby the possibility opens up to test the Cz/B 
and its interpretation as source of Neutrino Dark Energy. 

The paper is organized as follows. In Sec. 1 we discuss the MaVaN scenario and 
focus on a determination of the time dependence of the neutrino masses. Sec. 2 reviews 
the treatment of the absorption of an EHECz/ by the Cz/B in terms of the damping factor 
comprising the thermal background effects. Furthermore, in order to make contact to 
treatments in the literature [30,31,33,35-37], we include in the discussion common 
approximations [33,35,36] which neglect part or all of the dependence of the damping on 
the relic neutrino momenta. Moreover, we extend the complete analysis to incorporate a 
possible variation of the neutrino masses with time. In Sec. 3 we present and compare our 
results for the survival probabilities of EHECz/'s with varying and constant masses which 
encode the physical information on all possible annihilation processes on their way from 
their source to us, again taking into account the thermal motion of the relic neutrinos. 
In order to gain more physical insight, in addition, we disentangle the characteristic 
features of the absorption dips caused by the mass variation by switching off all thermal 
background effects. Sec. 4 illustrates the discovery potential of neutrino observatories 
for the Cz/B and gives an outlook for the testability of the MaVaN scenario. Therefore, 
both for astrophysical sources and for a topological defect scenario, we calculate the 
expected observable EHECz/ flux arriving at earth which results from folding the survival 
probabilities with the corresponding EHECz/ source emissivity distribution. In the latter 
case, we perform the full state-of-the-art calculation with the help of fragmentation 
functions and by the inclusion of all thermal background effects. In Sec. 5 we summarize 
our results and conclude. 

2. Mass Varying Neutrinos (MaVaNs) 

In Ref. [10] a new non-Standard Model interaction between neutrinos and a light 'dark' 
scalar field, the so-called acceleron, was introduced. In essence, it serves as possible 
origin of the apparent accelerated expansion of the universe and promotes the Cz/B 
to a natural dark energy candidate. Furthermore, as a very interesting and intriguing 
secondary effect, it causes a time evolution of neutrino masses. 

A follow up publication [15] takes care of a possible stability problem of the 
model [16-18] and furnishes a viable model of the whole scenario. 

Largely following Refs. [10,15], in this section we discuss the details of the complex 
interplay between the acceleron and the neutrinos that arises from a Yukawa coupling 
between them. Thereby, we will mainly focus on the determination of the resulting time 
variation of neutrino masses to be implemented later on in our analysis on relic neutrino 
absorption. For the latter it will turn out that the results are largely independent of 
the details of the model, since only a few generic features of the setting enter the 
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investigation. 

The new acceleron- neutrino interaction has a twofold effect. On the one hand, as a 
direct consequence, the neutrino masses m n are generated from the vacuum expectation 
value (VEV) A of the acceleron and become functions of A, m u .(A), i = 1,2,3. On 
the other hand, the dependence of m Vi on A turns the neutrino energy densities p Ui 
into implicit functions of A, since the energy densities p v . (m u .(A)) depend on the 
masses m Ui (A), i = 1, 2, 3. In this way, the energy density contained in a homogeneous 
background of neutrinos can stabilize the acceleron by contributing to its effective 
potential V e s(A). In other words, the dependence of the free energy on the value of 
A gets a contribution from the rest energy in neutrinos in addition to the pure scalar 
potential V (A). The total energy density of the system parameterized by the effective 
acceleron potential Vrff(vA) takes the following form, 



This is to be contrasted with the situation in empty space: if Vo(-4) is a 'run-away 
potential', the acceleron does not possess a stable vacuum state but rolls to its state of 
lowest energy given by the minimum of its pure potential Vo(.4). 

Taking now the expansion of the universe into account, the dilution of the neutrino 
energy densities p Vi [z) ~ (1 + z) 3 introduces a time dependence (here parameterized in 
terms of the cosmic redshift z) into the effective acceleron potential V e s- Consequently, 
in the adiabatic limit the equilibrium value A of the acceleron has to vary with time 
in order to instantaneously minimize its effective potential V^fr(^4)§. Finally, as the 
neutrino masses m u .(A) are directly affected by changes in the A condensate, they are 
promoted to dynamical quantities m Vi ( z ) depending on cosmic time. 

Note that Eq. ((2)) takes the neutrino energy density p Vi to be spatially constant. 
To justify this assumption, the A condensate is not allowed to vary significantly on 
distances of the order of the inter-neutrino spacing r of the relic neutrinos, with currently 
1/r ~ 336 1//3 cm _1 , where we have assumed a neutrino and anti- neutrino number density 
of n Uoi = n„ 0i ~ 56 cm -3 per species i = 1,2,3. Consequently, remembering that the 
range of the force mediated by a scalar field is equal to its inverse mass, one arrives at 
an upper bound on the A mass given by < 1/r ~ O(10~ 4 eV) at the present 
time. 

Let us now determine the time evolution of the physical neutrino masses m u .(z). 
Since the neutrino masses arise from the instantaneous equilibrium value A, we have 
to analyze the minimum of the total energy density V^(y4.). Assuming —^-^ to be 

| Under the assumption that the curvature scale of the potential is much larger than the expansion 
rate, d 2 V c s(A)/dA 2 — m 2 ^ H 2 , the adiabatic solution to the equations of motion applies. In this 
case for |^4| < Mpi ~ 3 x 10 18 GeV the effects of the kinetic energy terms can be safely ignored [10]. 
§ Since therefore in the presence of the relic neutrinos the acceleron possesses a stable (time dependent) 
vacuum state, in the literature both the acceleron and its VEV are referred to as A. 
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non-vanishing, one arrives at 



dV eS (A) ^dp Ui (m 



E 



dA 4f dm Vi 
i=i 1 



OA OA 

m v ,=m Vi {A) 



where [64] 



^ (1+ ^/w±3 and (4) 

vr 2 v 1 J ey + 1 v ' 

o 

t?v 

with = 1.69 x 1CT 4 eV denoting the neutrino temperature today. Note that the 
condition for the minimal energy density leads to a dependence of the neutrino masses 
on the neutrino energy densities which evolve with z on cosmological time scales. 

The smallness of the active neutrino masses m Ui can be explained by letting them 
only indirectly feel the acceleron via the seesaw mechanism [11-14]. Therefore, following 
Refs. [10,15], we introduce three 'right-handed' or 'sterile' neutrinos Ni with no Standard 
Model charges, whose masses Mjv 4 are constructed to vary with A due to a direct Yukawa 
interaction. In the seesaw mechanism the active neutrino masses m Ui are functions of 
the sterile neutrino masses M N .(A). Consequently, the A dependence of the M N .(A) is 
transmitted to the active neutrino masses m Vi {A) and causes them to change accordingly. 
Let us consider the interaction [10,15]: 

C D m Di .NiP h + KijNiNjA + h.c. + V {A). (6) 

where i,j = 1,2,3 are the family-number indices and correspond to the left-handed 
active neutrinos. Furthermore, kA corresponds to the A dependent mass matrix of 
the sterile neutrinos and tud is the Dirac type matrix (originating from electroweak 
symmetry breaking). Assuming the eigenvalues of kA to be much larger than the 
eigenvalues of mo one can integrate out the sterile neutrinos Ni , arriving at the following 
effective low energy Lagrangian [10, 15], 

£ D Mij(A) v h v h + h.c. + V (A), where (7) 

Mij(A) = (m ^ mD) ^ (8) 

represents the mass matrix of the active neutrinos. 

In order to solve Eq. (jnj) for m Vi (z) and to do MaVaN phenomenology the 
fundamental scalar potential Vo(-4) has to be specified in an appropriate way. Namely, 
the coupled neutrino acceleron fluid has to act as a form of dark energy which is stable 
against growth of inhomogeneities [16] and, as suggested by observations, must redshift 
with an equation of state u> ~ — 1 today. 

An appealing possibility arises in the framework of so-called hybrid models [65]. 
Those models were introduced to explain accelerated expansion in the context of 
inflation. In essence, two light scalar fields interact in such a way that one of them 
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stabilizes the other one in a metastable minimum. The energy density stored in the 
potential associated with the false minimum can drive accelerated expansion. 

It turns out, that a straightforward supersymmetrization [15] of the MaVaN model 
naturally sets the stage to apply the idea of the acceleration mechanism to dark energy. 
Identifying the acceleron with the former of the two light scalar fields, the hybrid model 
provides a microscopic origin for a quadratic acceleron potential V (^4) ~ A 2 . The 
role of the residual light scalar field coupled to the acceleron is attributed to the scalar 
partner Af of a sterile neutrino naturally present in a supersymmetric theory. The 
acceleron is stabilized by the presence of the fermionic neutrino background which 
drives its VEV to larger values. Accordingly, acceleration lasts as long as the VEV 
of the acceleron is sufficiently high to keep the sterile sneutrino Af in a false metastable 
minimum. So long as the energy density stored in the A condensate is sufficiently 
small, the combined scalar potential V(Af, A) will appear as dark energy redshifting 
with an equation of state u ~ — 1 [15]. Consequently, the neutrino dark energy density 
tlx ~ const, cosmologically behaves very much like a cosmo logical constant. 

According to Ref. [15] naturalness arguments require Af = Afi, assigning the lightest 
sterile sneutrino Af\ to be responsible for dark energy. Furthermore, one can conclude 
that Af\ has to be at least moderately relativistic today {m U01 < T uo ). Accordingly, 
in this supersymmetric MaVaN model probable instabilities [16-18] of highly non- 
relativistic MaVaN theories do not occur. 

In the past, the heavier two sterile sneutrinos A/2, A/3 of the theory were stabilized 
by the acceleron like the lightest Ai. However, by today they are assumed to have 
reached their state of lowest energy having acquired vacuum expectation values. || 

The relevant contribution [15] to the superpotential is given in terms of couplings 
of the superfield containing the acceleron A to two superfields with generation indices % 
and j, with i, j = 1, 2, 3, which comprise the sterile neutrinos N^Nj and their respective 
scalar partners Afi and Afj. The coupling constant matrix has elements Kij. This 
superfield interaction provides the necessary couplings mentioned above, namely of the 
scalar acceleron to the sterile sneutrino fields as well as the scalar acceleron Yukawa 
coupling to the sterile neutrinos in terms of k^. In Ref. [15] the one loop radiative 
corrections were estimated and taken to give a lower bound on the natural size for the 
magnitudes of the soft susy breaking masses squared of the A and the Afi- 

Given the exploratory nature of this investigation and following largely Ref. [15] 
it is reasonable to exploit the rough proportionality Sm 2 ^. ~ -m^ to get an estimate 
of the A mass squared m 2 A ] here Smjy. represent the one loop radiative corrections to 
the mass of a sterile sneutrino Afi and mn i label the respective eigenvalues of the Dirac 
type matrix mr>. For simplicity, we assume no off-diagonal elements for the coupling 
constant matrix k and denote the diagonal matrix elements by k^. Finally, one arrives 

I We refer to the mechanism proposed in [15] to evade large A/2,3 VEV contributions to the A mass 
conflicting with the upper mass bound set by the current inter-neutrino spacing O(10~ 4 ) eV. 
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(9) 



i=l 



such that the quadratic acceleron potential can be expressed in terms of neutrino mass 
parameters according to 

V,{A)^m\A 2 = j^KlmlA 2 . (10) 



i=i 



Now we are in a position to determine the respective mass-redshift relations m Ui (z) 
of the active MaVaNs whose mass squared differences today have to be compatible with 
oscillation experiments. Taking the acceleron Yukawa matrix k as well as the Dirac 
type mass matrix mo in Eq. (jHJ) to be diagonal, one arrives at the approximate seesaw 
formula for the physical neutrino masses: 



m Ui (A) 



ni 



Di 



K;A 



where i = 1,2, 3. 



Accordingly, the instantaneous minimum of V e R in Eq. Q is determined by 



dV cS (A) 
dA 



E 



dp Vi {m vv z) 



dm Ui 



_ m c, 



-m 



Di 



KjA 2 



2j2^rn 2 D A = 0. (12) 



i=l 



Since this equation has to hold for all z, the acceleron VEV becomes a function of z. 
As a direct consequence, it generates redshift dependent neutrino masses m Vi {z), 



m Vi {z) 



which implies 



m Vi {z) = m vo . 



KiA(z) 

A(0) 



with m^(0) = m^ Q 



A( z y 



(13) 



(14) 



Note that in general a MaVaN mass with subscript has to be identified with the present 
day neutrino mass. Consequently, the m Uo . have to be consistent with the mass squared 
differences measured in neutrino oscillation experiments. 

Accordingly, inserting Eq. (fTIJ) and Eq. (fTTH) into Eq. yields, 



E 

i=l 



m 



dp v .(m Ui ,z) 



rn 



m v .=m. 



Mm 

0,i A(z) 



0. 



m,, 



Evaluating Eq. (fTHjl at z = 0, 

3 6 3 

Z—/ rr)2 Z—/ 



m 



dp Vi (m u 0,i) 



i=l u o,i 



i=l 



dm. 



(15) 



(16) 



nv, =rrc„ n 



,0 



allows to eliminate 2 — from Eq. (|THjl leading to 

i=i m "o,i 



E 

i=i 



n? 



"0,-. 



dp v .{m Vi ,z) 



dm v . 



-4(0) 
"0,i .A(z) 



M(*)V ^(m^,0) 



0. (17) 



m v .=m. 



0,i 
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Finally, the solution for (^jy) > which can only be determined numerically fixes the 
neutrino mass evolution according to Eq. (|14|) in terms of the present day neutrino 

masses m UQi . 

However, the mass behavior in the low as well as in the high redshift regime can be 
approximated analytically by using the respective limits for J^p-. As mentioned before, 
in the supersymmetric MaVaN model [15] the lightest neutrino is at least moderately 
relativistic today such that its mass has to be very small, m Voi < T Vo = 1.69 x 10~ 4 eV. 
Furthermore, it can be deduced from the mass squared differences measured at neutrino 
oscillation experiments (see [66] for a recent review) that the heavier two neutrinos are 
non-relativistic today (m Uoi 3> T Vo for i = 2,3). 

Accordingly, in the low redshift regime it is a good approximation to neglect the 
contribution of the lightest neutrino species to V e s and solely employ the non-relativistic 
limit of J^- with i — 2, 3 and Xi = T ^> 1, 

oo 

9 ^^,z) ^ s f_t_ dy _ (18) 



dm Ui 7r 2 J e y + 1 

o 

Accordingly, in the low redshift regime Eq. (JT7j) is solved by 

A low (z) = A(0) (1 + z) (19) 

- m„. low (z) = m VQ .{l + z)-\ i=l,2,3, (20) 

where Eq. (|T4*Jl was used. 

Once in the past all neutrinos were relativistic. In this regime, Xj -C 1, J^- can be 
approximated by, 

oo 

dp^z) ^Tl (i + zf ^ [_!J dy (21) 



dm Vi 7T 2 v ' 1 J ev + 1 

o 

By taking the appropriate approximations Eq. (J 18)) and Eq. (|21|) for the two terms in 
Eq. ([T7]). V e s is minimized for, 

Aigh(^) oc (1 + zfl 2 (22) 
-> m Uthi Jz) oc (1 + z)- x l\ i=l,2,3, (23) 

where the factor of proportionality is a function of the present day neutrino masses and 
the integrals in Eq. (fTHj) and Eq. (j2Tj) . 

In our analysis on cosmic neutrino absorption later on, these approximations 
will help towards a better understanding of the numerical calculations since the 
corresponding results agree very well in the respective regimes. 

As becomes apparent from the approximations in Eq. (|20|) and Eq. ()23|) the MaVaN 
masses m Ui (z) are decreasing functions of redshift. Correspondingly, MaVaNs can be 
regarded as practically massless in the past whereas in the present epoch they have 
reached their maximal mass values m UQi . In the following we will assume a normal 
neutrino mass hierarchy and take the mass of the lightest neutrino to be 

m, 01 = 10- 5 eV, (24) 
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Figure 2. In this log-log plot the exact mass-redshift relations m Ui (z) are plotted as 
solid lines. In the low as well as in the high redshift regime they are well approximated 
by simple power laws, m Ui {z) oc (l+z) -1 andm^z) oc (l+z) -1 / 2 , respectively (dashed 
and dotted lines). The lightest MaVaN is assumed to have a mass of m„ 1 = 10 -5 eV 
today. Consequently, for a normal neutrino mass hierarchy, solar and atmospheric 
mass splittings fix the present time neutrino masses of the heavier neutrinos to be 
m„ 2 = 8.3 x 10" 3 eV and m Vo 3 = 5.17 x l(T 2 eV. 



such that it is still relativistic today (i.e. m VQ1 < T UQ ). According to the solar and 
atmospheric neutrino mass splittings the corresponding present time masses of the 
heavier non- relativistic neutrinos are fixed to be (see e.g. [66]) 

m V02 = 8.30 x lCT 3 eV > T vo , (25) 
m„ 03 = 5.17x KT 2 eV>TV (26) 

Figure El shows that the exact mass-redshift relations m Ui (z), i = 1,2,3, as numerically 
determined from Eq. (fTTj) in combination with Eq. (|14j). are very well approximated in 
the low as well as in the high redshift regime by simple power laws stated in Eq. (J2Uj) 
and in Eq. ()23j) . respectively. These results have to be contrasted with the standard 
time-independent neutrino masses m Uo . for i = 1,2, 3. 

Let us summarize the essential features of the presented viable MaVaN model which 
will enter the analysis on relic neutrino absorption later on. Firstly, the lightest neutrino 
is required to be still moderately relativistic today (m U01 < T Vo = 1.69 x 10~ 4 eV) and 
therefore fixes the neutrino mass scale to be low. Secondly, the specific mass-redshift 
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evolution m n (z) is determined by the model dependent quadratic form of the potential 
Vq ~ A 2 in Eq. (JTUJl which enters Eq. (|16|). The masses behave like oc (1 + z)~ x and 
oc (1 + z)~ x l 2 in the low and in the high redshift regime, respectively (cf. Eq. (j2Qj) and 
Eq. ()23|) ). However, the generically important feature (of any standard MaVaN model) 
is that the MaVaN masses m Ut (z) are decreasing functions of redshift which is easily 
achieved within the framework of the seesaw mechanism (see above). As will turn out 
in our investigation later on, for decreasing m Vi (z) the results do not strongly dependent 
on the actual functional dependence m Vi (z). Note that the characteristic differences in 
the absorption features of MaVaNs with respect to constant mass neutrinos (cf. Sec|U 
and SecEJ) are independent of the supersymmetrization of the MaVaN scenario. 

3. Damping of Extremely High-Energy Cosmic Neutrinos Traveling 
Through the Cosmic Neutrino Background 

Due to the feebleness of the weak interaction cosmic neutrinos can propagate 
cosmological distances through the cosmic microwave and neutrino background (CMB 
and Cz/B) without a significant probability of interacting. 

An interesting exception arises under the assumption that EHECz/'s with energies 
of order E$ = M 2 z /2m Vi = 4.2 x 10 14 GeV (O.OleV/m^) in the rest system of the 
target v exist, where Mz is the Z mass and m u . the respective neutrino mass. The 
Z resonance in the s channel for the process vv — > X, characterized by the energy 
E^f , enhances the cross-section for the annihilation of such an EHECz/ on a big-bang 
relic anti-neutrino (and vice versa) by orders of magnitude. As a consequence, the 
corresponding interaction probability significantly increases with respect to non-resonant 
scattering. Accordingly, the annihilation processes would reduce the survival probability 
of EHECz/'s traveling through the Cz/B to us and could therefore be detectable as 
absorption dips in the EHECz/ spectra. 

In this section we consider EHECz/'s (on whose plausible sources we will comment 
in Sec. EJ) propagating through a thermal bath of relic neutrinos in the expanding 
universe. Following Ref. [37], we discuss the corresponding damping rate of the EHECz/'s 
which governs their survival probability (cf. Sec. HJ. Furthermore, we will summarize 
common approximations for the damping which result from averaging over the neutrino 
momenta [36] or from completely neglecting the relic neutrino motion [30,31,33,35]. 

The investigation applies to both constant mass neutrinos as well as to MaVaNs. 
In the latter case, the neutrino masses m Vi are not constant but - as discussed in Sec. El 
- complicated functions of the neutrino densities and thus functions of z (cf. Eq. (JT7j) 
in combination with Eq. (|14j) in Sec. EJ). 

For simplicity, throughout this section we drop the indices i = 1,2,3 labeling the 
mass, energy or momentum of the neutrino mass eigenstates, since the discussion applies 
to all three neutrinos likewise (m u corresponds to m Ui , E to Ei etc.). 

The crucial quantity which describes the attenuation of an ultra-relativistic EHECz/ 
neutrino traversing the Cz/B is the damping rate j^. It can be expressed in the following 
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form^f, 



oc 



dp 



lvf> {E) = I — P 2 MP,T U ) a vi? (P,E), (27) 



and is governed by the Z-formation cross-section a u p(P, E) weighted by the momentum 
distribution f„(P, T u ) of the relic anti-neutrinos, both to be discussed in the following. 
Here, P denotes the modulus of the relic anti-neutrino momentum, whereas E represents 
the EHECz/ energy and T v is the CzaB temperature. Note that in the expanding universe 
these quantities, as well as the relic neutrino energy E p = \/P 2 + rrv*, are subject to 
cosmic redshift. They can be expressed in terms of their present day values, labeled by 
a subscript 0, in the following way 

P = P (l + z), E = E (l + z), T v = T VQ {1 + z) and E p = E po (l + z), (28) 

where we take the CzaB temperature at present to be T uo = 1.69 x 10~ 4 eV. 

Let us first consider the momentum distribution of the relic (anti-) neutrinos. Since 
they were ultra-relativistic at the time of decoupling their distribution assumes the 
relativistic Fermi-Dirac form [67] both for constant as well as varying neutrino masses. 
Neglecting the chemical potential of the relic (anti)-neutrinos it is given by 

U(P, T v ) = UP, T v ) = eP/ l + l - (29) 

Note that P/T v = Pa(l + z)/T„ d (l + z) is a non-redshifting quantity where the subscript 
d denotes the parameter values at decoupling. Accordingly, the relic (anti-)neutrinos 
can expand freely, preserving the form of f p (P,T u ) in Eq. (|2~9"j) . + 

In order to express the cross-section a u p(P,E) of the vv annihilation process it is 
convenient to introduce the parameter £ = Y\jM\ <C 1. It describes the square of the 
quotient of the total width for Z decaying to fermion pairs, Tz = 2.4952 GeV [68], and 
the mass of the Z, M z = 91.1876 GeV [68]. Accordingly, a uP (P,E) can be written in 
the following form [37], 

my ' ^2 2£J2 PE P J,_ ( S _M|) 2 + { S 2' ^ ' 

where Gp = 1.166 37 x 10 -5 GeV~ 2 is the Fermi coupling constant and E p = \JP 2 + ml 
is the energy of the relic neutrino. Furthermore, s is the square of the center-of-mass 
energy of the neutrino-anti-neutrino system. Using \J E 2 -\- m 2 v ~ E for an ultra- 
relativistic EHECz/ one arrives at the following expression for s, 

s = 2ml + 2E(E p - Pcos9), (31) 

The formula for the cross-section as well as the damping rate in Eq. (|27|l apply likewise if both of 
the interacting particles are replaced by their anti-particles. 

+ Strictly speaking, the neutrino distribution is not an equilibrium one, since more precisely the Fermi 
factor in Eq. © reads exp(^/(P/T iy ) 2 + [m v {T Vd )/T Vd ) 2 ) [64]. However, since T Vd - 1 MeV and 
accordingly m ^ d ^ < io~ 6 for m u (Td) < 1 eV, the bulk of the neutrino distribution with J^- > ""jj*^ 
is not affected by the strongly suppressed mass correction (neither in the case of constant mass neutrinos 
nor for MaVaNs) and can thus well be characterized by an equilibrium temperature T v according to 
Eq. ®. 
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where 6 characterizes the direction of the relic anti-neutrino with respect to the line of 
flight of the incident EHECz/ in the center-of-mass system. Accordingly, for fixed P and 
E the integration over s corresponds to the angular integration. As a consequence, the 
limits of the integral in Eq. (|3U|) take the values, 

s± = 2ml + 2E(E p ±P), (32) 

corresponding to cos 6* = =pl. 

Note that, following Ref. [37], in Eq. (J37?j) the energy dependence of the Z boson 
width r(s) in the Z propagator due to higher order corrections [69,70] has been taken 
into account. In the region of the resonance where the s dependence is crucial, T(s) can 
well be approximated by the linear relation [69,70] 

r s 

F ( s ) = T§2 S = VCiT' where Tz = const -- ( 33 ) 

±V±z IVlz 

The integral in Eq. (f3TJ|) for the cross-section can be solved analytically. Since the MaVaN 
mass m v {z) is a decreasing function of redshift, it has its maximal value today which 
corresponds to the mass of a constant mass neutrino. Accordingly, for both constant 
mass neutrinos and MaVaNs one can exploit that m v <C M z , E for all redshifts such 
that one finally gets the following expression for the cross-section [37]* 



c-uv{P,E) 



2V2G F T Z M Z 



+ 



+ 



2EE„ 



Ml 



1 



AEP(1 + 2 
1-C 



In 



1 + 04E 2 (E P + Pf - AM 2 z E{E p + P) 



Ml 



[1 + £) 2 AEPTz 



1 + 04E 2 (E P - Pf - 4M 2 E(E p -P) + M* 
f2E(l + 0(E p + P) 



arctan 



Mf 



T Z M Z 



arctan 



( 2E{l + Q{E p -P) 



M 2 



T z Mz 



(34) 



The calculation of the damping defined by Eq. (|27j) in combination with Eq. (|3*4*j) 
includes the full thermal background effects. To allow for a comparison of our findings 
to published results, in the following, we will summarize two common approximations 
used in the literature. They result from averaging over the neutrino momenta [36] or 
from completely neglecting the relic neutrino motion [30,31,33,35]. As will become 
apparent in the next section, in the neutrino mass range favored by the supersymmetric 
MaVaN model (cf. Sec. |2J) their applicability is quite limited. However, the assumption 
of relic neutrinos at rest will later on remove the thermal distortion of the absorption 

* The neglection of the neutrino mass in Eq. I|34|) leads to a spurious singularity in the integrand of 
Eq. (|27|l in a region of the relic neutrino momenta where the integral is supposed to be negligible. This 
can be cured by an appropriate limitation of the integration interval. 
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dips and thus allow a deeper insight into the characteristic features caused by the mass 
variation Tn u .(z) described in Sec. El 

The weaker approximation for a u p(P, E) commonly used in the literature is obtained 
by approximating the mean value theorem. The factor (s — 2m 2 u ) in the integrand 
of Eq. ()30j) takes the role of the weight function and is integrated over, whereas the 
residual part of the integrand is - as an approximation - evaluated at the midpoint of 
the integration interval s ~ 2EE p = 2E^JP 2 + m 2 . Accordingly, the cross-section takes 
the following form 

a u9 (E,P) = ^(s) = 2V2G F T z M z 8 (35) 

(s-Mf) +£s 2 

with £ = r|/M§. Let us point out again that due to higher order corrections a precise 
handling of the Z propagator near the resonance [69,70] leads to an energy dependence 
of the Z boson width T{s) (cf. Eq. (|33p). Note that in earlier work on neutrino- 
absorption (e.g. [30,31,33,35]) this correction has not been made but the lowest order 
(simple Breit-Wigner) form for the cross-section a uP (E,P) has been used. We would 
like to stress that by the averaging procedure, which leads to the approximation for o vD 
in Eq. (|35p. part of the angular information gets lost. This results in an underestimation 
of the thermal spread of the absorption dips [37]. 

Let us now come to the second, more radical approximation: In earlier work on 
cosmic neutrino absorption (e.g. [33,35]) it was assumed that the relic neutrinos are 
at rest [33,35], thereby switching off all thermal background effects. The corresponding 
cross-section can be recovered from the full expression in Eq. (J3*U|) by taking the limit 
P — > or from eq. (|33j) by setting P = such that s = So = 2Em v . In this case, 
the remaining integral over P in Eq. ()27|) reduces to the neutrino number density 
n v {z) = n p (z) = n u0 (l + z) 3 . Accordingly, the damping j v9 takes the following form 

-,UE) = nM = 2V2G F T Z M Z ^ - | ^-^ p - - ^ „„. (36) 

As will become apparent in the next section, by assuming the neutrinos to be at rest 
one neglects two conspiring effects on the damping which become more important with 
decreasing ratio m u /T u [37]: 

On the one hand, the full cross-section a vP (P,E) in Eq. (j2U), which governs the 
damping, depends on E p = P 2 + m 2 . As a consequence, the peak of the cross-section 
for a thermal bath of relic anti- neutrinos at E\ es [36], 

M 2 

E- s = . , (37) 

2(^P 2 + m 2 Ui -P t cosd) 

actually lies at lower energies than the one of a u p(so) for relic anti-neutrinos at rest: the 
energy E 0)i reduces to 

M 2 

Er = 5^. (38) 
Note that in the case of MaVaNs the masses m„, are functions of redshift, m v (z), and 



therefore introduce a z dependence into the resonance energies. Thus, they only coincide 
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with the respective constant mass ones for z = and m v (Q) = m VQi . We will discuss 
the consequences for the absorption features in detail in the next section. 

As indicated by Eq. ()H7j) . the effect of the relic neutrino momenta P becomes 
significant for small neutrino masses, according to Ref. [37] for m u < 0.01 eV. 
Furthermore, a„p(so) overestimates the peak height of the full expression and cannot 
account for the broadening of o~ u p(P, E) for increasing relic neutrino momentum P. 

On the other hand, the thermal distribution of the relic neutrinos which gives rise 
to a Fermi momentum smearing of the cross-section is totally neglected. In the full 
expression the damping results from the integration over all neutrino momenta, 
where the weight factor P 2 f p (P,T u ) selects relic neutrino momenta P of the order of 
T u . Accordingly, Eq. (}3l)j) overestimates the damping efficiency with respect to the full 
expression defined by Eq. (J34"|) and Eq. (J2JJ). As we will see in the next section the 
realistic description of neutrino-absorption leads to less well defined absorption features 
spread over a larger range of EHECz/ energies than in the idealized scenario which 
neglects any thermal effects. These discrepancies increase with decreasing neutrino 
mass. 

4. Survival Probability of Extremely High-Energy Cosmic Neutrinos 
Traversing the Cosmic Neutrino Background 

The relevant quantity to be discussed in this section is the survival probability P Va of 
extremely high-energy cosmic neutrinos v a of flavor a = e, fi,r traveling through the 
Cz/B to us. It is governed by the damping rate 7„p introduced in the last section and 
it determines, folded with the respective EHECz/ source emissivity distribution the 
resulting neutrino spectra to be observed on earth, which are treated in Sec. El 

The main goal of this section is to work out the characteristic differences in the shape 
of the absorption dips in the EHECz/ survival probabilities which arise from considering 
the neutrino masses to be dynamical quantities m Vi {z) instead of constant parameters. 
After presenting our results and pointing out the generic differences, we will have a 
closer look at the MaVaN absorption features. As motivated in the last section, for 
the purpose of gaining more physical insight, we will disentangle the different influences 
which define their shape. First of all, we will study the impact of the mass variation 
m Vi (z) as well as of the cosmic redshift caused by the expansion of the universe. To this 
end, we will initially switch off any thermal background effects by assuming the relic 
neutrinos to be at rest and compare the results for MaVaNs to those of constant mass 
neutrinos. 

Nonetheless, we would like to stress again that due to the low neutrino mass scale 
required in the MaVaN model under consideration (cf. Sec. EJ) only the full treatment 
of the background effects can serve as a test for Neutrino Dark Energy. 

In our calculation we make the standard simplifying assumption that the EHECz/ 
source switched on at a distinct redshift z s in the past. As concerns plausible EHECz/ 
sources, in the following we would like to mention the most relevant classes as well as 
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the corresponding typical EHECz/ emission redshifts z s . 

As a first possibility, EHECz/ are assumed to originate from pion decays, where 
the latter either are produced in inelastic pp or pry interactions. Those astrophysical 
acceleration sites (bottom- up mechanism), notably active galactic nuclei (AGN) and 
gamma-ray bursters, have source positions z s of a few (e.g. [38,71]). The conjectured 
energies of cosmic neutrinos produced by these astrophysical acceleration sites in the 
case of shock acceleration are E m3uX < 10 11 — 10 12 GeV [38,71]. However, even higher 
energies are possible in proposed non-shock acceleration mechanisms, such as unipolar 
induction, acceleration in strong magnetic waves in plasmas (wakefields) [72], or by 
magnetic recombination in the vicinity of massive black holes [73, 74] (see [71] for a 
recent review). 

As a second possibility, extremely energetic cosmic neutrinos with energies above 
10 12 GeV may be generated in the decomposition of so-called topological defects (top- 
down scenarios) into their constituent particles. Topological defects are predicted to 
originate from symmetry breaking phase transitions immediately after (hybrid) inflation 
(see [75] for a recent review). In particular, cosmic string formation is highly generic 
in Supersymmetric Grand Unified Theories (SUSY GUTs) (see [76] and references 
therein). Those topological defects produce super- heavy quanta generically denoted 
as X particles (often heavy Higgs or gauge bosons) with masses mx ~ 10 12 — 10 16 
GeV. Those X particles rapidly decay into stable Standard Model (SM) particles, also 
releasing extremely energetic neutrinos [39,77] with energies up to ~ 0.5 mx [78]. For 
those exotic, non-accelerator sources, z s can be as high as the epoch of light neutrino 
decoupling, z s ~ C(10 10 ) [32]. 

In our analysis we take resonant Z-production caused by the interaction with the 
Cz/B as the only source of attenuation of the propagating EHECz/. This approximation 
is well justified in the energy regions of the absorption dips which we are focusing on in 
our investigation [35] . Accordingly, the survival probability of an extremely high-energy 
cosmic neutrino z/j with i = 1, 2, 3 injected at redshift z s is given by, 

(/: -7^(£o(l + z)) 



P Vi (E 0} z s ) = exp 



1,2,3, (39) 



H(z)(l + z) 

where the integral in the exponential, which governs the survival probability, is called the 
optical depth (or the opacity). It contains the product of the propagation distance dr = 
dz/[{l + z)H(z)] and the damping rate ^ V y{E) defined in Eq. (|2*7jl with E = E (1 + z), 
which is integrated over all redshifts from the present time up to the emission redshift 
z s . Moreover, in a universe with negligible radiation component the evolving Hubble 
factor is given by 

H(z) = H (n M (i + zf + n k (i + zf + n A ) 1/2 . (40) 

As suggested by recent observations [68], we take a present day matter density VL M = 0.3, 
a curvature density Q k = and vacuum energy density = 0.7 as corresponds to a 
Lambda Cold Dark Matter (ACDM) universe. Note that this specific form for H (z) also 
applies to the MaVaN scenario under consideration: according to Sec. 121 the neutrino 
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Figure 3. Flavor survival probability P Va , a — e,fi,T with all thermal background 
effects included and integrated back to z s — 5 (left panel) and z s = 20 (right panel), 
respectively, for a normal neutrino mass hierarchy with m vo 1 = 10 -5 eV and varying 
neutrino masses (solid lines) as well as constant neutrino masses (dotted lines). 



dark energy density Qx redshifts with an equation of state lu ~ — 1 [15] and therefore 
behaves very much like a cosmological constant A, Qx ~ ~ const.. 

We will express our results for the survival probabilities in terms of the propagating 
neutrino flavors v a according to, 



in \'P vo with 



i 

1, 2, 3 and a = e, /i, r, 



(41) 



where the absolute square of the leptonic mixing matrix elements U a i relates the neutrino 
flavor components v a to the mass eigenstates z/«. Note that since the mixing matrix 
element \U e ?\ 1, the absorption dip produced by the heaviest mass eigenstate will 
not be visible in the case of P Ue . However, apart from this exception, the flavor survival 
probabilities to be discussed in the following exhibit absorption dips at the respective 
resonance energies of the mass eigenstates z/j, for i — 1,2, 3. 

Let us begin our analysis by comparing our results for the flavor survival 
probabilities defined in Eq. (|3*^j) and Eq. (|41j) . each with varying and constant mass, 
including all thermal background effects according to Eq. (|2Tj) and Eq. Throughout 
this section we assume a normal neutrino mass hierarchy, where the masses to take values 
according to Eq. (}2"2|) - Eq. (f2TI|) and the mass variation is given in Fig. |2 (cf. Sec. |2J). 



Probing Neutrino Dark Energy with Extremely High-Energy Cosmic Neutrinos 



20 



For numerical values of \U a i\ we take [79,80] (1 

/ 0.84 0.54 0.08 \ 
\U ai \ = 0.44 0.56 0.72 . (42) 
v 0.32 0.63 0.69 j 

Let us start by considering the flavor survival probabilities P Va for an emission 
redshift of z s = 5 which are plotted in Fig. El on the left hand side. Apparently, both in 
the case of varying and constant neutrino masses, the absorption features produced by 
the lighter two neutrino mass eigenstates are merged together. However, for MaVaNs, 
the absorption dips produced by the heavier two neutrino mass eigenstates exhibit sharp 
spikes at the respective resonance energies E^f = 2 ^ z with i = 2, 3 instead of being 
washed out and distorted to lower energies as in the constant mass case. In principle, 
by this means, the neutrino masses m„ i for i = 2,3 could be directly inferred from the 
respective spike positions of the MaVaNs dips in the low redshift regime. 

Altogether, the MaVaN absorption dips are much deeper and narrower in 
comparison to the constant mass features. In addition, the respective minimum positions 
are shifted to higher energies by almost an order of magnitude with respect to the 
corresponding constant mass dips. 

As demonstrated by Fig.|3]on the right hand side, for an increased emission redshift 
z s = 20, the absorption features are considerably deeper than for z s = 5. In addition, 
for each neutrino flavor the dips have merged, both in the case of varying and constant 
neutrino mass. As we learned in the last section, this effect of the thermal motion has 
increased with z, since the thermal bath of relic neutrinos was hotter at earlier times. 

For MaVaNs, the characteristic narrow spikes at the resonance energies E^f with 
i = 2,3 are less pronounced than for z s = 5 and also suffer a distortion towards 
lower energies. Nevertheless, for and v T , they remain well distinguishable from the 
respective constant mass dips. As in the case of z s = 5 the absorption features are 
clearly shifted to higher energies and exhibit substantially deeper dips. 

The characteristic absorption features produced by the mass variation can be 
worked out by separating the different influences on the MaVaN absorption dips. Let 
us for this purpose assume the relic neutrinos to be at rest, in order to eliminate any 
thermal background effects on the MaVaN survival probabilities. By doing so, we are 
left with the combined effects of the cosmic redshift and the mass variation. Let us first 
of all consider the former effect which is present both for MaVaNs as well as for constant 
mass neutrinos. It originates in the expansion of the universe and manifests itself in 
an energy loss of EHECz/'s of energy E according to E = E/(l + z), where E Q is the 
EHECz/ energy to be measured at earth. Accordingly, the survival probability P n of an 
EHECz/ is reduced, as long as somewhere on its way to us (z s > z > 0) it has the right 

jj Apparently, unitarity is not exactly fulfilled for the best fit values in Eq. Ij42[) . This leads to a small 
shift in the flavor survival probabilities P Va , which we have compensated for in FigEI~ FigElin such a 
way that P Ua — 1 outside the energy regions of the absorption dips produced by the mass eigenstates 
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Figure 4. Approximated flavor survival probability for P Va , a = e, [1,t, which assumes 
the target relic-neutrinos to be at rest, P = 0, after an integration back to z s = 5, for 
a normal neutrino mass hierarchy with m VQ 1 — 1CP 5 eV and varying neutrino masses 
(solid lines) as well as constant neutrino masses (dotted lines) plotted as a function of 
their energy E at earth. 



amount of energy, 

^res 

(ifk) " E ° £ (43) 

to annihilate resonantly with a relic anti-neutrino (whereas for all other energies the 
CVB is transparent for the EHECz^ such that P Va = 1). As a consequence, the effect of 
cosmic redshift is observable in a broadening of the EHECz/ absorption dips. 

Let us stress that in addition to this cosmological effect in the case of MaVaNs, the 
variation of the neutrino masses m u .(z) causes a redshift dependence of the respective 
resonance energies El es (z) as already mentioned in the last section. To be more concrete, 
the masses m Vi (z) at redshift z determine the corresponding resonance energies to be 
El cs (z) = M§ /2m Ui (z) in the rest system of the relic neutrinos which only coincide with 
Eft = Ml/2m„ . for z = 0. 

In Fig. |3] we plot the respective survival probabilities of the neutrino flavors with 
varying and constant masses, neglecting the relic neutrino momenta and integrating 
back to z s = 5. As expected from the discussion in the last section, the thermal spread 
of the absorption features provoked by the relic neutrino motion is removed such that 
the dips do not merge. Instead, for constant mass neutrinos, the absorption features in 
Fig. |U are only subject to the broadening caused by the cosmic redshift and span the 
energy interval specified in Eq. (j4Hjl . 
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In striking contrast, the absorption dips produced by the MaVaN mass eigenstates 

M 2 

exhibit sharp minima at the resonance energies E^f = 2m z with i = 1,2,3 showing no 
spread towards lower energies. Thus, interestingly, these MaVaN dips look like those of 
constant mass neutrinos in a non-expanding universe. Actually, it turns out that the 
mass induced redshift dependence of the resonance energies E\ es (z) = 2mv Z {z) ' * = ^> ^> ^' 
compensates for the energy loss of the EHECz/ due to cosmic redshift in the limit of low 
redshifts. This can be understood by remembering that the approximation in Eq. (J2(Jj) 
gives a good estimate for the redshift dependence of the neutrino masses m u . (z) in the 
low redshift regime (cf. Fig. EJ). In this limit one arrives at the following functional 
dependence of the resonance energies on z, 

M 2 

£ ™ = 2^T £ "« 1 + -~>- < 44 > 

In turn, the resonance energy El es (z) = + z) of an EHECz/ at redshift z 

corresponds to the redshifted energy measured at earth, 

EJ eB (z^ 



r| S = E™, with i = 2, 3, (45) 

Accordingly, the annihilation of an EHECz/ on the Cz/B at any given redshift z s > z > 
always leads to an absorption peak at E™?. In other words, in this approximation 
the square of the center-of-mass energy, so,i = 2m Ui (z)E (l + z) = 2m UQi E , i = 
1,2,3, becomes redshift independent. Correspondingly, the undistorted shape of the 
annihilation cross-section 0w(s Oj i) (weighted with the neutrino density per unit redshift 
and integrated over z) gets projected on the sky. 

Note that even for higher redshifts the appropriate mass-redshift approximation 
from Eq. (j2HJ) leads to the following redshift dependence of EJ eB (z), 

M 2 

E?*(z) = oc (1 + z) 1/2 with i = l,2, 3. (46) 

2m„ - " 



Apparently, the neutrino mass variation still partially counterbalances the effect of the 
cosmic redshift at high redshifts. 

In summary, all of the respective MaVaN absorption lines can be distinguished from 
those of constant mass neutrinos, since the redshift distortion is much less pronounced 
in any case. We would like to point out that this is a generic feature of any standard 
MaVaN scenario where the neutrino mass is a decreasing function of redshift. 

We are now in a position to complete the interpretation of the MaVaN features 
in Fig. El which result from combining all effects on the absorption features including 
the thermal ones. As already mentioned in the last section, the Fermi-weight factor 
P 2 f j? (P,T) entering the damping integral in Eq. (}2Tj) selects momenta P of the order 
of the Cz/B temperature T u (z). Accordingly, the ratio m v jT v (z) is a measure for the 
relevance of the mass and its probable evolution with respect to the temperature effects. 
Let us in the following discuss the absorption features which are completely determined 
by the thermal effects. In these cases the absorption lines were produced by relativistic 
neutrinos, since m Ut (z)/T u (z) 1. Firstly, for v e the absorption lines of MaVaNs and 
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constant mass neutrinos are similar (cf. Fig- EJ) - This can be understood by recalling 
that v e is mostly composed of the lightest mass eigenstate for which mass effects neither 
today nor in the past have played any role, since {m v%1 /T Vo <C 1). Secondly, as opposed 
to the case of z s = 5, the MaVaN absorption features for z s = 20 are distorted to much 
lower energies (cf. Fig. EJ). In addition, for MaVaNs, the low energy ends of the dips for 
and v T have the same shape as the one of v e . The reason is, that the absorption lines 
in this energy region stem from absorptions at high z < z s (as indicated by Eq. (}4"3*j) 
in absence of any thermal effects), where all neutrino masses still were negligible with 
respect to the temperature. As a consequence, they are clearly distinguishable from 
the corresponding ones of constant mass neutrinos, which apparently are already non- 
relativistic in the same energy region. This is due to the fact that while the temperature 
rises with increasing z, only the MaVaN masses m Vi (z) evolve and become lighter. Thus, 
MaVaNs generically turn non-relativistic much later than constant mass neutrinos. At 
energies above this transition from the non-relativistic to the relativistic regime, the 
variation of the heavier two neutrino masses is not washed out by the temperature 
effects. Therefore, it leads to sharp and thus deep absorption minima at the respective 
resonance energies according to Eq. (|45jl (cf. Fig. EJ), without and with transition of the 
regimes, respectively). 

As a conclusion we have learned that the characteristic effects of the neutrino mass 
variations in the case of the heavier two MaVaNs become apparent in the higher energy 
regions of the absorption dips, where the MaVaNs are still non-relativistic. However, also 
the low energy end of the absorption dips differs as long as the MaVaNs are relativistic 
and the constant mass neutrinos have already turned non-relativistic. 

The next section will deal with realistic neutrino fluxes to be measured by neutrino 
observatories. In general, a flavor tagging at extremely high energies seems unlikely 
apart from exceptions (e.g. for particular flavors), whereas all observatories will be 
sensitive to the flavor summed EHECz/ fluxes Yl F Va with a = e, [/,, r. Accordingly, we 
will base our final discussion on the totalized fluxes ^ F Va . In addition, we will include 
our results for the v T whose identification will at best be feasible by the LOFAR radio 
telescope [81]. For this purpose, for emission redshifts z s — 5, z s = 20 and z s = 50 
in Fig. El we collect both the resulting normalized sum of the survival probabilities, 
| Y Pui which governs the Y F Va as well as the respective P Ut . Note that emission 
redshifts of z s = 50 (and much higher) will contribute to the EHECV fluxes which result 
from the decomposition of super-heavy particles produced by topological defects (top- 
down sources). Apparently, also for emission redshifts of this order, the characteristic 
differences between the respective absorption features of varying and constant mass 
neutrinos, which we have discussed above, are still visible. Furthermore, the absorption 
dips get considerably deeper with increasing emission redshift z s . 
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Figure 5. The normalized sum of the survival probabilities, hY^^ni an d the flavor 
survival probability P Vt , respectively, including all thermal background effects, for a 
normal neutrino mass hierarchy with m„ 1 = 10 eV and varying neutrino masses 
(solid lines) as well as constant neutrino masses (dotted lines) plotted as a function of 
their energy E at earth for z s = 5, z s — 20 and z s =50. 



5. Absorption Dips in Realistic EHECz/ Spectra — Test of Neutrino Dark 
Energy? 

So far, in our analysis we have concentrated on the determination and the comparison 
of the survival probabilities P Va of mass varying and constant mass EHECz/'s traversing 
the Cz/B, where the P Va contain the physical information on the annihilation process 
v a v a — > Z. In the following we will outline how our results can be employed to perform 
relic neutrino absorption spectroscopy and to what extend the latter can serve as a test 
for the MaVaN scenario. For this purpose, we will firstly consider astrophysical EHECV 
acceleration sites and, secondly, topological defect sources both introduced in the last 
section. In the latter case, we will not approximate the EHECz/ injection spectrum by a 
standard power-law, but go through the appropriate calculation involving fragmentation 
functions as discussed below. 

As a starting point, we require both the detection of EHECz/ fluxes in the energy 
region of interest as well as the observation of absorption lines in these spectra. The 
EHECz/ flux F Ua (E ) for a neutrino of flavor a — e, //, r to be measured at earth is given 
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by [35] 




The flux integral is governed by the survival probability P v . defined in Eq. (f3*H|) weighted 
by the EHECz/ source emissivity distribution C U0 , which depends on the properties of 
the source as described in the following. On the one hand, the diffuse source emissivity 
£ Uf} (E (l + z),z s ) takes into account the distribution of the sources in the universe 
(the activity if) and on the other hand it considers the number of neutrinos of flavor 
(3 = e, /x, r emitted by each of the sources (the injection spectrum J v „). Under the 
standard assumption of identical injection spectra for all sources, one can factorize the 
z and E dependence, 

C„ (z s , E) = n(z s )J U0 (E), with E = E {1 + z s ). (48) 

As already stated in the last section, a flavor tagging at extremely high-energies 
cannot be expected at all neutrino observatories. However, we can hope for the 
identification of the v T absorption lines at LOFAR [81], which we will therefore consider 
according to Eq. (|47|). Furthermore, in our analysis, we will consider the total flux of 
all neutrino flavors F Ua which can well be approximated by [35] , 

oo 2 

J2 F ^( E o) - 4^ / ^lc^(E (l + z s ),z s )Y,P Vj (E Q (l + z s ),z s ), (49) 
o i =1 
where £* ot denotes the total, flavor-summed neutrino emissivity at the source and the 
formula holds as long as C Ufi + C Vt = 2£„ e . The latter is fulfilled for hadronic sources 
like astrophysical accelerator bottom-up sources or non-accelerator top-down sources, 
since in both cases the neutrinos emerge from charged pion decays such that 

C Ve : £,,. : I',,. 1:2: (I. (50) 

However, Eq. ()49|) also holds in the case of equal flavor source emissivities, 

C Ve : C M : C Vt = 1 : 1 : 1, (51) 

as could arise in the decays of topological defects not directly coupled to matter (e.g. 
mirror-matter 'necklaces') [82,83]. 

Note that the dependence on the leptonic mixing matric elements \U a j\, present in 
Eq. (JUJ), has dropped out in the expression for ~Y^F Va in Eq. (f4T?j) due to unitarity. 

In the next subsection we start our investigation by considering astrophysical 
(bottom-up) EHECz/ sources. In the subsequent subsection we continue our analysis 
for the case of topological defect (top-down) EHECz/ sources. 

5.1. Astrophysical neutrino sources 

In the following we will discuss EHECz/ fluxes which are assumed to originate from 
astrophysical EHECz/ sources. In order to parameterize their source emissivity 
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Figure 6. Projected sensitivity of LOFAR [48] expressed in terms of the diffuse 
neutrino flux per flavor, corresponding to one event per energy decade and indicated 
duration, together with E^F with F = J2^ a + E-^k, (l eit column) and EqF t with 
F T = F Vt + Fy T (right column) for varying (solid lines) and constant (dashed lines) 
neutrino masses and for z max = 20, assuming a normal neutrino mass hierarchy with 
m Vo l = 10~ 5 eV, n — 4 and a — 2 as well as £ max = 4 x 10 16 GeV. 

distribution C v „ we employ the following standard ansatz (e.g. [35, 84, 85]) in 
combination with Eq. ()48j) . 

V&s) = VoO- + z s) n 0{z s - z min )9{z m&x - z s ), (52) 
J U0 (E) = j U0 E- a e(E - E mhl )6(E m&x - E). (53) 

Throughout our analysis, we will take z min = and -Emi n = as default values and 
suppose that E max > E^(l + z max ) for i — 2, 3. Furthermore, we will not examine the 
possibility of broken power-law injection spectra, but assume the spectral index a to be 
constant in the whole energy region of interest. 

For the purpose of illustrating our results, we consider Eq. and Eq. for 
n = 4 and a = 2 as often used in the literature to mimic astrophysical sources and 
take E max = 4 x 10 16 GeV. As in the last sections, we assume a normal neutrino mass 
hierarchy according to Eq. (}2"4"|) - Eq. (}2l)j) and for the MaVaNs a neutrino mass variation 
according to Fig. El We present our results in Fig. EJ on the left hand side we plot the 
EHECz/ energy squared times the flavor summed flux EqF with F = J] Fv a + Yl ^k* 
according to Eq. ()49)1 and on the right hand side EqF t with F T = F Vr + F Pt as defined 
in Eq. (|47jl . In Fig. El we plot our results both for varying (solid lines) and constant 
(dotted lines) neutrino masses for z max = 20, together with the projected sensitivity of 
LOFAR [48] to be in operation by 2008 expressed in terms of diffuse fluxes per neutrino 
flavor, respectively. 

In our calculation we have assumed the EHEC^ flux to be close to the so-called 
cascade limit [86,87]. It applies to sources where the neutrinos emerge from pion 
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decays or even from electroweak jets [88] and are thus accompanied by photons and 
electrons which escape the source. Consequently, the measurements of diffuse gamma- 
ray fluxes, which are of the same origin but have cascaded down in energy during the 
propagation through the universe, restrict the neutrino flux to lie below the cascade 
limit. Apparently, the predicted sensitivity of LOFAR [48], corresponding to one event 
per neutrino flavor per energy decade, lies below the cascade limit by several orders 
of magnitude. Accordingly, at best 3500 neutrinos (plus anti-neutrinos) in the energy 
interval 10 12 — 10 13 GeV can be expected to be detected by the radio telescope. 

The large event numbers N result in tiny error bars (a = V3N /3iV). In a blow- 
up of the absorption features in Fig. we have adjusted them to the curves with no 
absorption for emission redshifts z max = 20, z max = 10 and z max = 5 from top to bottom. 
Thereby, we have assumed a conservative and therefore rather poor energy resolution 
corresponding to one energy bin per energy decade, whereas at best LOFAR is predicted 
to achieve an energy resolution of AE/E ~ 30% [81]. The latter would correspond to 
~ 4 energy bins per energy decade. Apparently, the dips become considerably deeper 
with increasing z max . Despite the underlying low neutrino mass scale, both for varying 
and constant neutrino masses, LOFAR can be expected to produce significant evidence 
for absorption dips in the EHECz/ spectra for emission redshifts 2 max = 20 and z max = 10 
- even for a bad energy resolution. In the case of 2 max = 5, in the interval 10 14 — 10 15 
GeV, the considerably higher dip depth for MaVaNs leads to a more than 5a deviation 
from the curve with no absorption whereas for constant mass neutrinos the departure is 
not significant (both for the flavor summed flux F and for F T ). Even if the underlying 
EHECz/ fluxes are much lower, at least for EHECz/ sources at z max = 20, a detection of 
absorption features produced by varying, light neutrino masses could well be feasible. 
Correspondingly, if such EHECz/ fluxes of astrophysical origin exist, the most direct 
detection of the Cz/B so far seems to be in reach within the next decade. From the 
experimental point of view, the prospects are even better for scenarios with time varying 
neutrino masses, which in general produce deeper absorption dips in the regime of 
astrophysical emission redshifts. 

Let us now turn in more detail to the prospects of probing scenarios of Neutrino 
Dark Energy by identifying the characteristic absorption signatures of a possible 
neutrino mass evolution. Since the EHECz/ fluxes are governed by the respective survival 
probabilities discussed in the last section, the characteristic differences in the absorption 
features for varying instead of constant neutrino masses are maintained. Namely, for 
MaVaNs one observes a clear shift of the dips to higher energies as well as considerably 
deeper absorption minima with respect to constant mass neutrinos. Accordingly, given 
a decent energy resolution of AE/E ~ 30% [81] for LOFAR, relic neutrino absorption 
spectroscopy could serve as a test for the nature of neutrino masses and therefore for 
Neutrino Dark Energy. However, certainly, the feasibility strongly depends on the energy 
resolution achieved by the EHECz/ observatory. 
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Figure 7. The energy squared times the flavor summed neutrino flux EqF with 
F = E Fv« + E f p q (left column) and E$F T with F T = F„ T + F„ T (right column) 
for varying (solid lines) and constant (dotted lines) neutrino masses for z max = 20, 
z max = 10 and z m ax = 5 from top to bottom, respectively. All curves assume a 
normal neutrino mass hierarchy with m vo 1 = 10~ 5 eV, n — 4 and a = 2 as well as 
£max = 4 x 10 16 GeV. 
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5.2. Topological defect neutrino sources 

In the following we will discuss neutrino fluxes expected to result from exotic top- 
down EHECz/ sources like topological defects. As already mentioned, EHECz/'s might 
be produced among other Standard Model particles in the decays of super-heavy X 
quanta which constitute the topological defects. Accordingly, the corresponding EHECz/ 
injection spectra J in Eq. ()48|) are fragmentation functions which can reliably be 
predicted by the help of Monte Carlo generators [89, 90] or via the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) evolution [91-93] from experimentally determined 
initial distributions at a scale Mz to the ones at mj. The corresponding injection rate 
(the activity i] in Eq. (J48J1). which in particular determines the overall normalization of 
the neutrino flux, in terms of cosmic time t is given by, 



where Qq is the energy emitted per unit volume per unit time at present and p is 
a dimensionless constant. Both Q and p depend on the specific topological defect 
scenario [77]. In our analysis, we calculate the absorption features in the EHECz/ 
spectra for topological defect sources by the help of fragmentation functions as well 
as by Eq. (j^lj) according to Eq. (J4*%j) .ff 

In the following, we will illustrate the prospects for performing relic neutrino 
absorption spectroscopy for top-down sources considering as example superconducting 
strings whose evolution is characterized by p = in Eq. ()54|). Note that in the energy 
region of the absorption dips the Standard Model and SUSY fragmentation functions 
(which we have taken from [92]) have practically the same shape. Accordingly, all our 
results on relic neutrino absorption are independent of the supersymmetrization of the 
MaVaN scenario (cf. Sec.^J). 

In Fig. |H1 we plot the expected absorption features for varying (dashed lines) and 
constant (dotted lines) neutrino masses as defined in Eq. (|2~4"j) -Eq. (|2fij) . where for the 
MaVaNs we assume a mass variation according to FigJ21 Again, we present our results 
for the energy squared times the flavor summed flux EqF according to Eq. ()49)1 in the 
left column of Fig. |H] and EqF t as defined in Eq. (}4"TJ) in the right column of Fig. |HJ As 
in the case of astrophysical sources, in the first panel of each column we have plotted 
the projected sensitivities of ANITA and LOFAR as well as the predicted EHECz/ flux 
for mx = 10 16 GeV with and without absorption (where the latter by design scratches 
the cascade limit). In the second panel of Fig. EJ we again show a blow-up of the 
absorption features. Apparently, with respect to the astrophysical sources (cf. Fig. EJ), 
the dips are broader by almost an order of magnitude in energy. This is due to the 
fact that the constituents of topological defects have started to decay and therefore to 
release neutrinos at z 1. As a further consequence, the dips for top-down sources are 

ffln the literature on absorption dips so far the injection spectra of top-down sources had been 
approximated by Eq. |(52Jl and Eq. J53J). Cf. e.g. Refs. [35,37] for the appropriate values for n 
and a. 




(54) 
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Figure 8. The energy squared times the neutrino flux EqF with F = J] Fu a + S ^v a 
(left column) and EqF t with F T = F Vt + F„ T (right column) for varying (dashed lines) 
and constant (dotted lines) neutrino masses expected from the decomposition of a 
superconducting string with p = and mx — 10 16 GeV, in the first panel together 
with the projected sensitivities for ANITA [42] and LOFAR [48], which correspond to 
one event per flavor, energy decade and indicated duration, respectively. 



much deeper than for bottom-up sources, since the survival probability of a neutrino 
traveling to us is much lower (as discussed in the last section and as illustrated by 
the comparison of the flavor summed survival probability for z s = 50 and z s = 20 in 
Fig. El). Consequently, both of these features facilitate a detection of the absorption dips 
in the EHECz^ spectra. Accordingly, top-down sources with the same underlying fluxes 
as astrophysical sources are even better suited to provide evidence for the existence of 
the CuB. By this means valuable information on the topological defect scenario, on 
cosmological parameters as well as the neutrino mass scale could be gained both for 
varying and constant neutrino masses. Furthermore, for the first time in cosmic particle 
physics, the GUT energy scale 0{mx) ~ 10 16 GeV could be probed. 



Probing Neutrino Dark Energy with Extremely High-Energy Cosmic Neutrinos 31 

In the second panel of Fig. |Sl we have also included the expected error bars, again 
assuming one energy bin per energy decade, both for the MaVaN and constant mass 
neutrino absorption lines. Clearly, the discrepancy between the two curves is larger for 
EqF t than for EqF, whereby in contrast to EHECz/s of astrophysical origin, constant 
neutrino masses produce somewhat deeper dips than time dependent masses. These 
features can be understood by realizing that v T is dominantly composed of the heaviest 
mass eigenstate m U3 , whereas EqF by definition gets equal contributions from all mass 
eigenstates. Furthermore, for constant mass neutrinos, m Vz is the only mass eigenstate 
for which the ratio m U3 /T u (z) ~ (1 + z) 3> 1 up to z ~ 1000. In other words, it produces 
much deeper absorption dips than the lighter mass eigenstates (even when integrating 
back to z ^> 1) and their characteristic shape is not dominated by the temperature 
effects (cf. the discussion in the last section). In contrast, for MaVaNs, the ratio 
m Vi {z)jT v {z) with i = 1,2,3 for all mass eigenstates drops much faster with increasing 
z and takes values m Ui (z) /T u (z) <C 1 for z 3> 1. As a result, the low energy end of the 
dip (which corresponds to higher annihilation redshifts z) has the same shape both for 
all MaVaN mass eigenstates and the lighter two constant ones and is totally determined 
by the thermal background effects. 

In summary, promisingly, a resolution of absorption features for either mass 
behavior seems to be possible both for ANITA and for LOFAR. Yet, a differentiation of 
the MaVaN and constant mass neutrino absorption features seems only feasible, if tau 
flavor tagging and a good energy resolution are achieved. 

6. Summary and Conclusions 

In light of the number of extremely high-energy neutrino (EHEC^) observatories in 
operation and under construction with a combined sensitivity ranging up to 10 17 GeV, 
the prospects for establishing the existence of EHECz/ fluxes appear to be very promising. 
As a next step, the exciting possibility opens up to trace the annihilation of EHECz/'s 
and relic anti-neutrinos (and vice versa) into Z bosons by localizing absorption dips 
in the EHECz/ spectra at energies set by the neutrino masses. On the one hand, their 
detection could furnish the most direct evidence for the Cz/B so far and thereby confirm 
standard cosmology back to the time of light neutrino decoupling. On the other hand, 
the shape of the absorption lines could reveal a variation of neutrino masses with time 
and thus verify the interpretation of the CiaB as source of Neutrino Dark Energy. 

We therefore considered a viable Mass Varying Neutrino (MaVaN) model with the 
following features entering our analysis on relic neutrino absorption. By the requirement 
that the lightest neutrino still has to be moderately relativistic today the neutrino 
mass scale is set to be low, which leads to very conservative predictions. Furthermore, 
the evolving neutrino masses m Vi (z), which we determined numerically as functions of 
redshift assuming m VQ1 = 10~ 5 eV, turned out to be well approximated by simple power 
laws (1 + z)~ l and (1 + z) _1//2 in the low and in the high redshift regime, respectively. 
Accordingly, as a generically important feature, they are decreasing functions of redshift 
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as in all standard MaVaN scenarios. 

In order to provide all technical tools to interpret EHECz/ absorption dips for a given 
injection spectrum and to extract valuable information on neutrino physics, cosmology 
and possibly physics beyond the Standard Model, we proceeded in the following way. 
We considered in parallel the neutrino masses to be functions of cosmic time as well as 
to be constants. In our analysis we took into account the full thermal background effects 
which result from the relic neutrino motion according to their phase space distribution. 
In order to compare our results to the literature, we included in our discussion common 
approximations [33,35,36] which neglect part or all of the dependence of the damping 
on the relic neutrino momenta. 

On the level of the survival probabilities which govern the EHECz/ fluxes, we 
found the following results: For low emission redshifts (z ~ 0(5)), the absorption dips 
produced by the varying neutrino masses m Ui (z) for i — 2,3 exhibit narrow absorption 
minima, which do not suffer a distortion to lower energies as the corresponding dips 
of constant mass neutrinos. As a consequence, for MaVaNs, the absorption dips of 
the flavor components and v T (which are mostly composed of the heavier two mass 
eigenstates) are clearly deeper and shifted to higher energies by almost an order of 
magnitude with respect to the corresponding constant mass minima. For an increased 
emission redshift z ^> 5, these features become somewhat less pronounced but essentially 
prevail. A better understanding of the characteristic signatures caused by the mass 
evolution was obtained by switching off the superposing thermal wash-out caused by 
the relic neutrino motion. After neglecting the relic neutrino momenta for this purpose, 
we found that the crucial deviations result from the dependence of the corresponding 
resonance energies on the neutrino masses E\ es ~ 1/m^ for % = 1,2,3. In the case of 
MaVaNs, the mass variation m u . (z) induces a dependence on the annihilation redshift z 
according to E] cs (z) ~ E^(l + z) for all neutrino species % = 1, 2, 3 in the low redshift 
regime. Accordingly, this z dependence of the resonance energies compensates for the 
energy loss of the EHECz/ due to cosmic redshift proportional to (1 + z)~ l , resulting 
in narrow absorption spikes at constant energies El es (z)/(1 + z) — E^ (like one would 
expect for constant neutrino masses in a non-expanding universe). In contrast, for 
constant neutrino masses m Vi = m UOi the absorption dips are broadened, since the 
redshifted resonance energies to be measured on earth are given by El es /(1 + z) — 
Eq^/(1 + z), for z taking values between and the EHECz/ emission redshift. Since in the 
standard MaVaN scenario the neutrino masses are decreasing functions of redshift, they 
generically reduce the effect of cosmic redshift on the EHECz/ survival probabilities. As 
a result, they always produce deeper absorption minima, which, in addition, are shifted 
to higher energies in comparison to the dips caused by constant neutrino masses. 

In order to illustrate the discovery potential for absorption dips in the EHECz/ 
spectra to be observed at earth and to estimate the prospects of testing scenarios 
of Neutrino Dark Energy, we considered plausible EHECz/ fluxes originating from 
astrophysical acceleration sites or from topological defect sources. We presented our 
results both for the energy squared times the flavor summed EHECz/ flux E%F with 
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F = J2F Ua + F Va and for EqF t with F T = F Ut + F 9t , where the latter can at 
best be identified by LOFAR [81]. Despite the adopted low neutrino mass scale, we 
found both for varying and constant neutrino masses that for topological defect and for 
astrophysical EHECz/ sources at £ max > 5, LOFAR and ANITA promise a statistically 
significant evidence for absorption dips (even if the underlying fluxes are well below the 
cascade limit). Accordingly, the most direct detection of the Cz/B so far seems to be in 
reach within the next decade. 

Furthermore, the flux dips of varying and constant mass EHECz/'s expected 
from astrophysical sources retain the characteristic differences induced by the survival 
probabilities. Besides being clearly shifted to higher energies, the MaVaN dips are 
deeper and therefore even facilitate a resolution of absorption features in the EHECz/ 
spectra in comparison to constant mass neutrinos both in the case of EqF and of EqF t . 
As a main result of our analysis, these deviations of the MaVaN and constant mass 
absorption curves for astrophysical sources turned out to be statistically significant, yet a 
decent energy resolution seems necessary for their detection. Given an energy resolution 
of AE/E ~ 30% as at best achievable for LOFAR [81], relic neutrino absorption 
spectroscopy could reveal a variation of neutrino masses and therefore possibly the 
nature of Dark Energy. 

As concerns topological defect sources, the absorption lines in the EHECz/ fluxes 
for time dependent and constant neutrino masses altogether are more similar in shape, 
however, somewhat deeper for constant neutrino masses. Furthermore, they extend 
to much lower energies than for astrophysical EHECz/ sources and their minima are 
considerably deeper. All of these features are a result of the much higher annihilation 
redshifts z s ^> 1 possible for EHECz/'s originating from the decomposition of topological 
defects in comparison to EHECz/'s from astrophysical acceleration sites. At high 
redshifts, the EHECz/'s are absorbed by a hotter bath of relic neutrinos. Consequently, 
in the energy region spanned by the absorption dips where m Vi /T v <C 1, thermal 
background effects wash out any features produced by the neutrino mass or its possible 
variation. Since the MaVaN masses are decreasing functions of redshift, they reach 
this limit for much smaller redshifts than the corresponding constant masses. Only the 
mass of the heaviest constant mass eigenstate is sufficiently large, m U3 /T u (z) 3> 1, in the 
relevant energy region, leading to a deeper absorption curve than the one produced by all 
of the other MaVaN and constant mass eigenstates. Since v T is mostly composed of the 
heaviest mass eigenstate, F T exhibits deeper constant mass dips than F. Accordingly, 
for F T the signatures of varying neutrino masses can more easily be distinguished from 
those of constant masses than in the case of F. However, in order to reveal a neutrino 
mass variation, it seems necessary both to identify the tau neutrino flavor and to have 
a good energy resolution. 

Recently, the authors of Ref. [17] claimed that certain constraints on the 
fundamental scalar acceleron potential Vo(*4) (cf. Sec.|2J can also lead to MaVaN models 
stable against the growth of inhomogeneities [16] even in the highly non-relativistic 
regime. Accordingly, one could construct other MaVaN models than the one under 
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consideration in this paper, whose viability would not rely on a low neutrino mass 
scale (cf. Sec. |2J). As we pointed out, the characteristic absorption signatures of any 
standard MaVaN scenario (cf. Sec. |U and Sec. EJ) are essentially generic apart from 
details. Yet, a higher neutrino mass scale would even increase the overall dip depth 
in comparison to our rather conservative predictions and also reduce the importance 
of the thermal background effects on the absorption features. Accordingly, we would 
expect the features produced by a possible mass evolution not to be washed out by 
the temperature effects for a wider energy range of the dips. Thus the deviations with 
respect to the corresponding constant mass curves would even be more prominent. 
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